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Abstract 

The following lectures are an introduction to the phenomena of par- 
tonic saturation and nonlinear evolution equations in Quantum Chromo- 
dynamics. After a short introduction to the linear evolution, the problems 
of unitarity bound and parton saturation are discussed. The nonlinear 
Balitsky-Kovchegov evolution equation in the high energy limit is intro- 
duced, and the progress towards the understanding of the properties of 
its solution is reviewed. We discuss the concepts of the saturation scale, 
geometrical scaling and the lack of the infrared diffusion. Finally, we give 
a brief summary of current theoretical developments which extend beyond 
the Balitsky-Kovchegov equation. 



1 Introduction 

One of the most intriguing problems of Quantum Chromodynamics is the growth 
of cross sections for hadronic interactions at high energies. Let us consider the 

tPresented at XLIV Cracow School of Theoretical Physics, June 2004, Zakopane, Poland. 
'Permanent address. 
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Figure 1: Scattering of two hadronic probes at high energy. Left: QED-type 
diagrams; right: diagrams with gluon self-interactions. 

scattering of two particles at very high energy as shown in FigQ As the energy 
grows, so does the probability of emission of soft particles. In the case of QED 
one has to consider a diagram of the type shown on the left hand side graph in 
Fig-E In QCD one also encounters these diagrams, but there are also additional 
diagrams of the type shown on the right hand side of Fig. ^ Since gluons are 
the carriers of the color charge and couple to each other, the increase of energy 
causes a fast growth of the gluon density and, consequently, of the cross section. 
This increase leads to the formation of a dense, colored medium at very high 
energies. 

In perturbative QCD the growth of the gluon density in the limit of high 
energy is governed by the BFKL Pomeron equation IT]. The solution to this 
equation gives a very strong, power-like growth of the gluon density, and of the 
resulting cross section 

fix) ^ x-^ , 

where x is the Bjorken variable (a fraction of target's longitudinal momentum 
) and A = '"^^'^ q,^ ig the intercept of the perturbative Pomeron in the leading 
logarithmic (LLx, Inl/x ^ 1) approximation. In the pioneering paper Gri- 
bov. Levin and Ryskin pointed out that the gluon recombination is important 
at high energies. It reduces the growth of the parton density by producing an ef- 
fect called the perturbative partonic saturation. In 2 , a new nonlinear evolution 
equation in double leading logarithmic approximation (DLL A, In l/xlnQ^ 3> 1) 
for the gluon density has been postulated. In addition to the gluon production 
it also takes into account the recombination effects 

^ dlnl/xdQ^ ' TT 3CFi?2Q2i^M^'y )J • (1) 

Note the negative sign in front of the nonlinear term responsible for the gluon 
recombination. The strong growth generated by the linear term is damped for 
large gluon densities xG{x, Q^) (of the order l/a^). 
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Partonic saturation plays also an important role in the context of the unitarity 
bound. It is well known that the hadronic cross sections obey the Froissart 
bound [3] which stems from the general assumptions of the analyticity and 
unitarity of the scattering amplitude. The Froissart bound implies that the total 
cross section does not grow faster than the logarithm squared of the energy 

CTtot = Ar {Insf , (2) 
ml 

where is the scale of the range of the strong force. It is generally believed, 
that the parton saturation mechanism leads to the unitarization of the cross 
section at high energies. Unfortunately, the problem is quite complex since the 
parton saturation is purely perturbative mechanism while the Froissart bound 
has been derived from general principles and it refers to the QCD as to a com- 
plete theory of strong interactions ( including the nonperturbative effects) . The 
GLR postulate resulted in increased efforts to develop a theory able to de- 
scribe the saturation at high energies. One effective theory for high density 
partonic systems at small x is the Color Glass Condensate |1| with the result- 
ing JIMWLK evolution equations Another approach has been developed 
by Balitsky [7] who constructed an infinite hierarchy of coupled equations for 
correlators of Wilson lines. In the mean field approximation the first equation 
of this theory decouples, and is equivalent to the Kovchegov equation :8, (de- 
rived independently in the dipole approach jSj). The Kovchegov equation is 
a nonlinear equation for the dipole scattering amplitude valid in the leading 
logl/a; approximation. The Balitsky-Kovchegov (BK) equation is perhaps the 
best known equation that includes the saturation effects and has a virtue that 
it can be relatively simply solved, at least numerically. 

The goal of the following lecture is to introduce the reader into the phe- 
nomenon of the partonic saturation and nonlinear evolution using BK equation 
as an example. We shall start with a brief review of the linear evolution in 
QCD. Then we shall recall the Froissart bound and the necessary conditions for 
its derivation. Then the properties of the solution of the BK equation will be 
investigated with particular emphasis on the the infrared diffusion, saturation 
scale and the geometrical scaling. We will continue with an analysis of this 
equation in general case in 4 dimensions, which takes into account the dipole 
spatial distribution in the impact parameter space. We will conclude with a 
short outlook and the discussion of recent theoretical developments in the field. 



2 DIS kinematics and variables 

Let us concentrate on the deep inelastic scattering process of Icpton off the 
hadron or nucleus. For completeness, let us first recall the basic kinematic 
features of DIS, as represented in Fig. El The total energy squared of the 
electron-nucleon system is S(^eN) = + whereas that of the photon-nucleon 
system is Sj*n = (p + q)'^- The photon virtuality is denoted hy — {k — k')^ = 
—Q^ < and the Biorken variable x = ^— — „i .'^ . 
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Figure 2: Deep inelastic scattering process of electron on the hadronic target. 
* * 





Figure 3: Left: photon with virtuality probes quark with a longitudinal 
momentum fraction x. Right: the virtuality (resolution) is increased so the 
density of quarks also grows. 



oo 



The high energy regime, is defined as 

S(7-Ar) 
X 

Y 



In 1 /x — > oo . 



0, 



(3) 



3 The linear evolution equations of QCD 

Let us consider a scattering of photon with virtuality off a hadron at center 
of mass energy y/s. The photon virtuality defines a resolution scale A ~ — 

with which one probes the partonic structure of a hadron, see the left hand side 
graph in Fig. |21 At a given resolution t = InQ^/Qg, the photon probes the 
density of partons q{x, t) with a fraction of the hadron momentum x = Q'^ /s. 
By increasing one also increases the resolution, and the density of quarks is 
larger: q{x, t) + 5q{x, t), see right-hand graph in Fig.|2| The underlying process 
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Figure 4: Additional set of diagrams present in the DGLAP evolution. 



can be described by the following linear evolution equation for density 
dq{x,t) as{t) /"^ dy 



dt 



y 



(4) 



The splitting function Pqq{z — x/y) describes the probability of finding a quark 
inside theparent quark, with a fraction z = x/y oi the parent quark momentum. 
This is one of the set of the well known DGLAP evolution equations ^Ol- For 
the DGLAP equations to be complete, apart from the quark density g(a;, t), one 
has to include the gluon density g{x, t) coupled to q. The evolution of the gluon 
density is shown in Fig. ^ The full set of singlet DGLAP equations reads as 
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where I](a;, i) — '}2,j\qi{x, t) + qi{x, t)] is the singlet quark density. 

An alternative approach to DGLAP X. is to consider a fixed virtuality of 
the probe and to increase the energy s (or alternatively rapidity Y), see Fig.|S| 
This procedure leads to the BFKL equation 

dG{x,t) asN, 



dt'JC{t,t')G{x,t') 



(5) 



d In 1/x TT 

which is an evolution equation in Bjorken x. The quantity IC(t,t') ( Lipatov 
kernel) describes the probability of branching of the gluon with virtuality t' into 
another gluon with of virtuality t. The function G{x, t) is called the unintegrated 
gluon density and is related to g{x,t) by 

g{x,t) = / dt'G{x,t') . 



Both /C and Py have perturbative expansions in as and share a finite number 
of common terms in the expansion. The solution to the BFKL equation has the 
form, see for example 



G{x) 



Xp = 41n2iVcQ;s/7r. 



(6) 



This solution strongly grows with 1/x (and, correspondingly with energy s 
Q'^/x) which contradicts the Froissart bound [3]. 
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Figure 5: Evolution in the BFKL framework: virtualities of probes are fixed 
and energy is increased. 

4 Froissart bound 

Let us recall the basic assumptions used in the derivation of the Froissart bound 
13 • This bound applies for the total cross section for scattering of two hadrons, 
and reads as follows 



CTTOT < ^(Ins) 



TO, 



(7) 



Its derivation is based on two assumptions: the unitarity of partial amplitudes 
and the finite length of the strong interaction. The first condition demands that 
the S matrix has to be unitary 

S^S = SS'' = 1. (8) 

A set of particle states |to) has to satisfy completeness relation 

^|TO)(m| = 1. (9) 



Then the probability that the initial state \i) evolves into the final state |/) is 

Pf^^\(f\S\^)\' . (10) 

Since a probability that a given final state comes from any initial state is 1 
therefore the sum of these probabilities over all the initial states equals 1 



J2Pf. = mS^S\z) = 1, 
/ 

which is interpreted as a unitarity condition for the S-matrix 



(11) 



The second assumption is on the finite range R of strong force determined 
by the mass scale to^ 

R^ — , (12) 

TO^ 

which determines the cutoff. The scale m-^ is entirely nonperturbative. The 
Froissart bound can be derived P] by using these two assumptions and the 
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Mandelstam representation. One has to stress that the Froissart bound is ap- 
phcable to the complete QCD theory that includes both perturbative and non- 
perturbative parts. 

It is worthwhile to mention that while one believes that Froissart bound is 
valid, it is not immediately visible from the data. The data show that the struc- 
ture function for deep inelastic scattering at high photon virtualities strongly 
rises with decreasing x, consistent with a power behaviour x~'^ with A ~ 0.35. 
Further, the data show no sign of a logarithmic dependence is seen, see Fig. O 




Figure 6: F2 structure function data from HERA collider and fixed target ex- 
periments. 

Pinning down the saturation effects in experimental data is a nontrivial task. 
A specific problem with the structure function data is that it involves completely 
inclusive measurements. In particular, structure function F2 is averaged over 
the impact parameter of the j* — p collision. It is also known that saturation 
effects crucially depend on this variable. Therefore the search for the saturation 
should involve more exclusive processes, e.g. diffraction. For more information 
on the saturation phenomenology consult | 12l \1'6\ and references therein. 

5 Parton saturation and nonlinear evolution 

Given that the Froissart bound should be satisfied, a natural question arises: 
how to modify the perturbative evolution in order to tame the growth of the 
cross section? One would like to identify the Feynman diagrams which are 
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responsible for gluon recombination and to derive the appropriate evolution 
equations that include these diagrams. As already stated in the introduction, 
the standard evolution equations lead to a strong rise of the parton density in 
high energy limits. One can expect that partons overlap when their density 
becomes large. The schematic picture of the saturation phenomenon is shown 
in Fig.0 






Figure 7: Schematic view of parton saturation. Horizontal axis is energy squared 
s, the vertical axis is r, the parton size. 



The horizontal axis represents the energy whereas and the vertical one the 
parton size defined by the inverse of the photon virtuality r ~ l/Q in the deep 
inelastic scattering process. The onset of saturation depends on energy and 
on the parton size. The larger the size of the partons, the earlier they fill up 
the available area and start to reinteract. A further increase of the energy 
will not increase the probed parton density. Therefore, apart from production 
diagrams, one has to include additional diagrams which take into account gluon 
recombination. This leads to the modification of the evolution equation by a 
term which is nonlinear in density. The first equation of this type was the GLR 
equation |2| 



^2 d^xGjx.Q') 
' dlnl/xdQ'^ 



asNc 



xG{x,Q^ 



3CfR^ 



2m2 



[xG{x,Q')] 



(13) 



The first term in H13|) is the usual DGLAP term in the double logarithmic 
approximation. In l/cc In Q^/A^ 3> 1, whereas the second, nonlinear, term is 
responsible for gluon recombination. The nonlinear term is inversely propor- 
tional to the hadron area ~ B? and the scale ~ at which the gluon density 
is probed. The smaller the hadron area, the earlier the partons fill it up and 
saturate. The scale defines the parton size r ^ For small values of r 

the saturation is delayed to larger energies. 

The GLR equation sums a subset of diagrams within the double leading 
logarithmic approximation. They are called fan diagrams, and are illustrated 
on the right hand graph in Fig. |H1 
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Figure 8: Left: linear evolution. Right: fan diagrams summed by the nonlinear 
GLR equation 

The nonlinear Balitsky-Kovchegov (BK) equation is valid in the leading log- 
arithmic In 1/x approximation. It has been derived independently by Kovchegov 
[5] within the dipole formulation of high energy scattering and by Balitsky \J\ 
from the operator product expansion for high energy scattering. More precisely, 
Balitsky's equations form an infinite hierarchy of coupled equations for correla- 
tors of Wilson lines, and only in the mean field approximation the first equation 
decouples and is equivalent to the equation derived by Kovchegov. An inde- 
pendent approach is that of Color Glass Condensate @] in which the evolution 
is governed by the JIMWLK functional equation "S* equivalent to Balitsky hi- 
erarchy. In this lecture we will study the solution of the Balitsky-Kovchegov 
equation which is currently the simplest tool to describe the parton saturation 
phenomenon. 

6 Multiple scattering in dipole picture 

In this section we will follow the derivation of the BFKL 9 and BK 'S' equations 
in the dipole picture. Consider a heavy quark- antiquark pair, onium, shown in 
Fig. El who's wave function in the momentum space is denoted by 



of light cone momentum carried by the quark. 

The dipole picture is formulated by going to the transverse coordinate space 




V'a/g (xo,xi,zi) 



/ 



C/3(ki,2i) , 



(2^)2 



e' 
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1-Zi p-kj 



Figure 9: Heavy quark-antiquark dipole onium. 

$(o)(xo,Xi,zi) = ^|^i'3(xo,xi,zi)|2 , 

where xq, xi denote the positions of the quark and antiquark respectively, which 
form the end points of the dipole. 

Then, as shown in Fig. ^| one adds one soft gluon with its fongitudinal 
momentum much smaller than that of the original quark (anti-quark) 22/^1 ^ 1 




l-z,,p-kj l-zi,p-kj 



Figure 10: Onium with an additional single soft gluon. 

The relation between the one-gluon wave function <i>*^^^ and onium wave 
function without any soft gluons <i>*^°^ is 

7^ Jzo ^2 J X^oX^2 

In the limit of large number of colors, the gluon can be represented by a quark- 
antiquark pair, as in Fig. 1111 The emission of one additional gluon is equivalent 
to the splitting of the original dipole (0, 1) into two dipoles (0, 2) and (2, 1) with 
probability of branching given by the measure 

X20X12 
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1-Zi ,P-k 1 

Figure 11: Oiiium wave function which consists of two dipoles. 



The process of emissions of subsequent soft dipoles can be repeated in the 
analogous way to obtain the wave function with an arbitrary number of gluons 
see Fig. El To describe such a process, Mueller |^ introduced a dipole 
generating functional 

Z(boi,Xoi,2:i,u) , 

which satisfies the normalization condition 

Z(boi,xoi, = 1) = 1 . 

The wave functions for arbitrary number of gluons can be obtained by perform- 
ing functional differentiation of Z, 

5 6 



a>(")(xo, 



Xl,X2, 



$(") 



-Z{xo,xi,zi,u) 



6u{x2) Su{x3) (5w(x„+i) 

Here ^^^"^ gives probability of finding n daughter dipoles that originate from 
parent quark-antiquark dipole (0,1). The daughter dipoles are produced in 
positions Xk with k = 2, . . . ,n. In the following, wc will use another notation 
with 

Xqi = Xo - Xl , 

representing the transverse size of the dipole and 



boi = 



Xo + Xl 



the impact parameter (position) of this dipole. 

By investigating the relation between wave functions with n and n+l dipoles 
Mueller derived [HI the following differential equation for the generating func- 
tional 

(iZ(boi,Xoi,?/,M) 



dy 



^20^12 



Z{hoi + ^,X20,?;, u) 2'(boi - ^,xi2,y,'u) - Z(boi, xqi, u) 



(14) 
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l-Zi,P-ki 



Figure 12: Oiiium wave function with arbitrary number of dipoles. 



where the evolution variable is the rapidity y — In 1/2:+. Using equation (|14|) . 
one can obtain the evolution equation for the scattering amplitude of dipole on 
the target. First, one constructs the dipole number density 

ni(a;oi,x,b- bo,y) = — -r- — -Z(boi, Xqi, F, , 
ou(h, x) 

and in general, the density for k dipoles 

nfc^nti ^ , zu^i. (15) 

The amplitude for scattering of single dipole on a target, see left graph in 
Fig. 1131 can be then obtained by convoluting the dipole number density with 
the propagator of that dipole in the nucleus 

A^i(xoi,boi,r) - J d[7'i]ni7i , (16) 

where d[V]i = ^:^(Phi is the phase space measure, and 7 = 7(x, b) is the 
propagator of a single dipole in the nucleus. By differentiating the equation for 
the generating functional and using the relation H16() one can obtain the linear 
evolution equation for the dipole-target amplitude 



dA^i(boi,xoi,y) 
dY 



rf^X2 xgi 



^20 ^12 



iVi (boi + ^ , X20 , y ) + A^i (boi - ^ , xi2 , y ) 



7Vi(boi,xoi,r) 



(17) 



It has to be stressed that only the contribution from the single scattering of one 
dipole on the target has been included in the derivation. The equation (|17|l is 
a dipole version of the BFKL equation in the transverse coordinate space as 
derived in W- One can also generalize this equation by taking into account a 
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multiple scattering of many dipoles on the target, see right graph in Fig.^l To 
this aim one takes the number density of k dipoles, Eq. 1)1511 . and then convolutes 
it with k propagators for respective dipoles. The following expression for the 
amplitude is then 



oo „ 

^(xoi , boi , F) = ^ / dlVk] Tifc 7i • ■ ■ 7fc 
k=i 



(18) 



where the measure is defined as 



By differentiation of the equation for the generating functional Z one can obtain 
the evolution equation for the amplitude which takes into account also multiple 
scatterings [S] 



dA^(boi,xoi,y) 
dY 



(P^2 Xqi 



^20 ^12 



iV(boi + ^, X20, r)+^(boi-^, xi2, r) 



7V(boi, xoi, F) - 7V(boi + X20, Y) JV(boi - Xi2, F) 



(19) 



The characteristic feature of this equation is its nonlinearity. Thus, in the 
dipole approach, the multiple scattering of many dipoles in the onium leads to 
nonlinear evolution equation for the amplitude. This has to be contrasted with 
the single scattering of one dipole which leads to the linear BFKL-type equation. 
One has to stress that this multiple scattering is a completely incoherent process: 




Figure 13: Left: single dipole scattering which leads to linear BFKL evolution 
equation (|17|l . Right: multiple dipole scattering which results in nonlinear 
Balitsky-Kovchegov evolution equation H19|l . 

dipoles scatter independently of each other and there are no correlations. This is 
quite an important simplification which results in a relatively simple and closed 
evolution equation. These correlations are now a subject of intensive research. 
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Figure 14: Schematic representation of the dipole position in impact parameter 
space. (xo,Xi) denote end points of the dipole. 



7 Balitsky-Kovchegov equation at high energies 

In the remaining sections of the paper will be devoted to the discussion of the 
solutions to the BK equation 



(i7V(boi,xoi,F) _ f d^xaxgi 



dY 



^20 ^12 



^^(boi + ^, X20, r)+A^(boi-^, xi2, r) 



7V(boi, xoi, r) - ^(boi + X20, Y) 7V(boi - xi2, r) 



(20) 



Let us highlit the salient features of this equation: 

i) BK equation is an evolution equation in rapidity y = In l/x. 

ii) It requires the initial conditions A^'^*'-'(boi, xqi, y = 0) which depend on the 
target of the specific process to be specified. 

iii) The BK equation is valid in the leading logarithmic approximation in which 
the powers in (a^ In l/x)" are counted. 

iv) In this approximation the strong coupling as is fixed. 

v) In (|20|l boi represents the impact parameter whereas Xqi is the size of the 
dipole, see Fig. ^1 The problem involves (4 + 1) variables: four degrees of 
freedom per dipole and one evolution variable. 



7.1 Toy model in (0 + 1) dimensions 

Searching for the solutions to Eq. H20() one has first to notice that this equation 
has two fixed points 

div(boi,xoi,y) _ 

dY 

which are at 

iV = and iV = 1 . 

It is quite instructive to first investigate the toy model in (0 + 1) dimensions, 
when amplitude depends only on the rapidity N = N{Y) and the kernel is 
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simply a constant. Then the equation reduces to 

^ = cj{N~N'),u;>0. 

The above equation was first discussed by Vcrhulst in 1838 as a model for self- 
limiting population growth in biology. The solution to this equation, called 
logistic curve, can be easily found 

^(^) - +Vi ' ^^"Hy = 0) = C, 

and is illustrated in Fig. 1151 Its crucial property is that it saturates for very 
large values of Y 

Vc^o N{Y) ""-^ 1 , 
in contrast to the solution of the linear equation which grows exponentially. 
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Figure 15: Illustration of the solution to the Verhulst equation (saturated line) 
and the linear equation (exponentially increasing). 

The toy model teaches us that the fixed point at = is unstable with 
respect to the linear part of the evolution. In contradistinction, the value of 
N = 1 turns out to be a fixed point. After sufficiently long evolution in Y the 
solution will reach this fixed point, independently of initial conditions (provided 
that iV(o)(y = 0)^0). 

8 Solution in (1 + 1) dimensions 

Having briefly looked at the toy model let us proceed to the full equation which 
depends on (4+1) variables which makes it very difficult to solve even numeri- 




15 




Figure 16: Rapidity dependence of the solution to (1 + 1) dimensional Balitsky- 
Kovchegov equation (saturated lines) for two, fixed values of the dipole size r as 
compared to the solution of the linear BFKL equation (exponentially increas- 
ing). 



cally. The biggest complication is that the variables boi and Xgi are entangled 
in the arguments of the functions iV, see Eg . 02011 . However, the kernel depends 
only on the sizes Xoi,X2o,xi2. By assuming that the solution N is translation- 
ally invariant 

7V(boi,xoi,r)^7V(|xoi|,r) , 

the problem is reducible to (1-1-1) dimensions with no dependence on the impact 
parameter boi. Physically, this corresponds to scattering of a dipole on an 
infinite and uniform nucleus. In this approximation, the BK equation in (1 -|- 
1) dimensions has been extensively studied numerically, 1151 1161 [T7| and 
analytically [HITniOni- 

In Fig. ^] we illustrate the rapidity dependence of N{Y,r — |xoi|) for two 
fixed values of r, and compare it to the solution of the linear equation. We 
notice that the solution to the BK equation has the same qualitative features 
as the toy model. For any given r the solution of the nonlinear equation tends 
to unity, whereas the linear solution increases exponentially. 

Also when the dipole size is larger, the system also saturates earlier, for 
smaller values of Y. To understand this feature better let us look at the solution 
as a function of the dipole size r for given, fixed values of Y. We notice that for 
smaller dipoles the amplitude = 1 as the rapidity increases. In the following, 
see Fig. 1171 we will also study the dependence of N on different initial conditions. 

In Fig. 1171 several initial conditions have been shown. They require differ- 
ent normalisations and exhibit different type of behaviour for large values of r 
{N{r) — > 1 or N{r) ^ as r — oo). In all cases the common feature of the 
solution is that TV = is an unstable fixed point and = 1 is the stable one. 
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Figure 17: Dipole size dependence of the solution to the (1 + 1) dimensional 
Balitsky-Kovchegov equation for different values of rapidity. From top to bot- 
tom: different initial conditions; left: linear scale; right: logarithmic scale. 
Dashed line denotes the initial distribution a,t Y — 0. Solid lines from right to 
left are for increasing values of rapidity. 
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8.1 Saturation scale 



The solutions shown in Fig. El exhibit three different behaviours : in one re- 
gion where the ampUtude is small the nonhnear corrections are negligible, the 
transition region and the asymptotic region where the amplitude N 1. The 
boundary of the transition region is characterized by a saturation scale Qs{Y) 



r < 



r > 



siY) 



QsiY) 



N totheleft , 
1 totheright . 



The saturation scale can be extracted from the solution to the BK equation. 
The leading rapidity dependence is exponential 

QsiY) = Qoexp(a,Ar)y-'^, A ~ 2.4 , 

with some subleading corrections |21[I19| . 

The qualitative properties of the BK equation are roughly similar to the 
properties of the dipole cross section of the Golec-Biernat and Wusthoff satura- 
tion model 13^ in which the following form of the dipole cross section has been 
postulated 



a(y,r) 



J (fhN{h,r,Y) = ctq 



1 — exp 



(21) 



In the saturation scale: Q^iY) — q0.28{y-Yo) rapidity dependent. The 
normalisation ctq been adjusted to obtain the best description of the experi- 
mentally measured proton structure function F2. In the regime where dipoles are 
smaller than the reciprocal of the characteristic saturation scale r < l/Qs{Y), 
the cross section is small and proportional to 

a{r,Y)/ao^r^Q^jY)/4, 

in accordance with the requirement of color transparency. When we consider 
large dipoles, r > l/Qs{Y), the cross section saturates to (Jq 

a(r,y)/ao~l , 

and becomes independent of both r and Y. 



8.2 Geometrical scaling and travelling waves 

While investigating the solution as a function of the dipole size r one observes 
that the solution reaches the universal shape independently of the particulars of 
the initial condition. For different values of Y the solutions have similar shapes 
but are shifted shifted towards smaller values of the dipole size. This property 
is known as geometrical scaling and it was first searched for and found in the 
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Q =: 



Figure 18: Schematic representation of dilute and saturated regions regions in 
the kinematic space {Q,Y). 



data at HERA electron- proton collider |22[ I23|. Mathematically, geometrical 
scaling means that the solution to the BK equation depends only on a single 
combined variable 

rQs{Y) , 

instead of r and Y separately, i.e. 

Nir,Y) = N{rQ,{Y)) . 

In terms of logarithms of variables, using the rapidity dependence of the satu- 
ration scale Qs{Y) ~ Qq exp{asXsY) one finds that 



Ini 



lng,(r) =lnr-f a.A.y 



If we interpret In r as a spatial coordinate and Y as the time variable, then the 
property of geometrical scaling implies that the solution is a wave front that 
propagates with a constant velocity asXs, see JHl- It has been also described 
as a soliton wave The scaling property is also present in the Golec-Biernat 
and Wusthoff saturation model (|21|l . 

The transition between the dilute and saturated regimes can be illustrated as 
in Fig.^l The dense and dilute regions are divided by the critical line identified 
as the saturation scale. The higher the rapidity the denser the system becomes 
and finally partons begin to reinteract. The saturation occurs earlier also, for 
the larger size of the partons. 



9 Diffusion properties of the BK equation 

Insofar we have looked at the solutions in the coordinate space in accordance 
with the original formulation of the BK theory. By Fourier transform to the 
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Figure 19: Left: solution to the BFKL equation in the momentum space as 
a function of momentum k for various fixed values of rapidity Y — 1, . . . , 10; 
right: the same but both BFKL (dashed) and BK (solid) solutions are showed. 



momentum space 

dr 



0(fc,y) := / —Joikr)N{r,Y) 







r 



in (1 + 1) dimensions one can obtain a more compact form of this equation 

^ a,J ^/C(fc, k') 0(fc', Y) - as<P\k, Y) . (22) 

In Eq. H22|l the integral operator JC{k,k') is the usual BFKL kernel in momentum 
space. 

The solution to the linear part is well known. In the saddle point approxi- 
mation it reads 



kcbik,Y) = \ exp(a,x(0)y) exp (~ !!'-^^y.^^l ) ■ (23) 

The first exponential is responsible for the fast increase of the gluon density 
with rapidity. The value of x(0) = 4 In 2 is the famous BFKL intercept. The 
second exponential causes the diffusion of the momenta into the ultraviolet 
and infrared regions. It is well known fact that the BFKL equation exhibits 
strong diffusion, which can be interpreted as a random walk in the In k space 
of transverse momenta. The rapidity (energy) plays here the role of the time 
variable. The left plot of Fig . 1 1 91 illustrates the numerical solution to the BFKL 
as a function of transverse momentum for fixed values of Y . 

The Gaussian shape predicted by H23|) is obvious and the width increases 
with rapidity. This consitutes a potential problem since while one starts from 
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a perturbative calculation at a fixed, large scale fco, eventually the nonpertur- 
bative regime of Aqcd ~ fc ^ fco is reached. This is not the case for the 
solutions nonlinear BK equation, which are illustrated by solid lines on right 
plot in Fig. ^1 Clearly, in the case of the solutions to the BK equation, the 
diffusion into the infrared region is strongly suppressed. With increasing Y the 
distribution peak moves away from the initial value fco towards the larger values 
of fc. Therefore one can define saturation scale as a position of the peak 

Q,(r) = fc„,ax(n . 

The suppression of the diffusion can be also visualized with help of the normal- 
ized distribution 

vi/(fc Y) = (24) 

On left hand side of Fig. [201 shows the contour plot in (fc, Y) space of ths dis- 
tribution H24() for the case of the linear BFKL equation. The contour lines 
correspond to the constant values of the normalized distribution \l/(fc,y). The 
diffusive character of the solution is clearly visible. On the right hand side of 
Fig. 1201 we show the corresponding contour plot for the nonlinear BK equation. 
We see that the contour lines are shifted towards the higher values of transverse 
momenta^. We can also identify a line in (fc, Y) space which divides the region 
where there still is a diffusion (to the right) and where there is no diffusion. 
The contour lines are parallel to each other in the latter case what means that 
the solution is scaling there. When the parametrisation ^ = Infc/fco — XY + ^o 
is introduced the solution depends only on ^ alone. The critical line defines 
the saturation scale that was introduced in the previous paragraphs. It turns 
out that the nonlinear BK equation can be approximately treated as a linear 
diffusion equation with the absorptive boundary close to the critical line defined 
by the saturation scale QsiY) This approximation allows us to evaluate 

precisely the rapidity dependence of the saturation scale. 

Recently, there has been quite substantial development achieved towards the 
understanding of the solutions to the (1+1) dimensional BK equation. In a series 
of important papers |19' it was proved that BK equation can be approximated 
as a diffusion equation with a nonlinear term. This makes it equivalent to the 
Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation [21] 

dtu{t,x) = dl{t,x)+u{t,x)[l-u{t,x)] , (25) 

where the change of variables from (Y, In fc) to {t, x) has been performed with 
simultaneous identification of (/> — s- u. FKPP equation has been already studied 
in many fields of physics, and its solutions are very well understood. In partic- 
ular it is well known that the FKPP equation has a travelling wave solution for 
large times (which are equivalent to large energies) which is just a property of 
geometrical scaling. For a review see for example 1251 . 

distortion of the contours at the highest values of k is unphysical and is caused by 
cutoffs in numerical calculation. 
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ln(k/k„) 



Figure 20: Contour plots of the renormalized distribution 'i>(k,Y) in the case 
of the hnear BFKL solution and the nonlinear BK solution. 

9.1 BK equation with running coupling 

As we have already stated, the BK equation has been derived within the leading 
logarithmic in In 1 /x (LLx) approximation in which the coupling constant is 
fixed. However, it is very well known fact, that NLLx effects are very important 
in the BFKL formalism At NLLx order the coupling runs, and the linear 
BFKL equation becomes very unstable. The reason is that the linear evolution is 
very sensitive to the details of the running coupling regularization. This effect is 
illustrated in Fig.l21lwhere we show the solution k(j)^^^^{k, Y) with the running 
coupling, as a function of the transverse momentum for increasing values of Y . 
For small values of rapidity the position of the maximum remains close to the 
value of the initial condition k ^ kg. However as the rapidity increases the 
position of the maximum abruptly shifts to the small values of momenta. The 
actual shape of the solution are critically dependent on the regularized value of 
the coupling, k ~ k^cg ^ fco- 

With the inclusion of the running coupling the BK equation has the following 
form 



The solution to the above equation is illustrated in Fig. superimposed 
onto the solution of the linear BFKL equation. The solution is much more stable 
than in the linear case, and the maximum of the distribution is not shifted to 
the infrared. Instead as rapidity increases, the maximum moves towards the 
higher values of transverse momenta. The reason for this is that the nonlinear 
term strongly damps the diffusion into the infrared regime. The saturation scale 
Qs{Y) provides a natural cutoff for the low momenta and then no dependence 



d(j){k,Y) 
dY 
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Figure 21: Solution to the BFKL equation with running couphng. 




Figure 22: Sohition to the BK equation (as compared to BFKL) with running 
couphng . 
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on the regularisation of the running couphng is seen (see for example |17[ 12 71 

EHlEni)- 

One might ask whether the geometrical scaling is still preserved in the pres- 
ence of the additional scale Aqcd which is implicitly introduced by use of the 
running coupling constant. It turns out that the scaling still holds, although 
the saturation scale has a different rapidity dependence (see ^| |^ 1301 1211 E] 
and also |2]) 



QsiY) = Aqcd exp (^^^{Y- Yo) + In^ Qo/AqcdJ , (27) 

with c ~ 2. The above formula has been derived by assuming that the local 
exponent of the saturation scale 

_ dln(Q,(y)/A) 
^^^^ dY 

has the similar form as in the case of the fixed coupling 

\{Y) - casiQliY)) . 

By using these two formulae one can derive the saturation scale for the running 
coupling l|?7jl . see pTHHT)) . 



10 Solution in (3 + 1) dimensions 

10.1 Spatial distribution: impact parameter dependence 

The phenomenon of saturation discussed so far has been based on the properties 
of the solutions to the BK equation in 1 + 1 dimensions. Such approach ignores 
the spatial distribution of the probe-target system. The solution to (1 + 1) 
dimensional BK equation is integrated over the impact parameter. This is a 
result of the fact that in 1 -I- 1 dimensional case we have assumed an infinite size 
of the target, and ignored any edge effects. One might expect that the more 
realistic picture of saturation looks as follows: as the probe collides with the 
target, a dense system of partons emerges in the limited region of the impact 
parameter space. For larger impact parameters, the density of partons becomes 
more and more dilute. The radius of the dense, saturated system expands with 
the growth of the energy. This process is schematically drawn in Fig. One 
expects that the impact parameter profile of the scattering amplitude has the 
behaviour shown in Fig. 1241 

A question thus arises whether BK equation can provide information about 
the impact parameter profile of the amplitude, and whether this profile is con- 
sistent with the qualitative picture of saturation. Let us recall the BK equation 
with full dependence on all coordinates 
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^ Total size of system 

Total size of system ^ f 



Size of dense system 




Figure 23: Schematic picture of the impact parameter dependence of the parton 
density at high energies. Right plot: as energy increases the area of the saturated 
region, the black disc, increases. 




b 



Figure 24: Impact parameter dependence of the scattering amphtude for two 
different values of rapidity. As energy (rapidity) increases the radius of the 
saturated plateau also grows. 
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Figure 25: Parametrisation of dipole position. 



djV(boi,xoi,y) 

dY 



I 





x: 



■20 -^12 L 



iV(boi,xoi,r) - 7V(boi + ^,x2o,r)A^(boi 




As stated before, in general the problem is very difficult, even numerically, 
since one has four degrees of freedom per dipole plus rapidity as the evolution 
variable. Even though the integral measure does not depend on the positions 
of the dipoles, the impact parameter dependence still exists. It is implicitly 
generated via the couplings of bij to xtj in the arguments of the functions iV, 
see Eq. 1281 To simplify the problem, and yet retain the information about the 
impact parameter dependence, we note that measure in the equation 



see Fig. 

Therefore one can assume that the position of the dipole is specified by 
three variables : the dipole size r, the impact parameter 6, and the relative 
orientation of the dipole with respect to the impact parameter axis (angle 9), 
see Fig. |2H1 The invariance with respect to global rotations is equivalent to a 
condition that the target is cylindrically symmetrical. Thus the problem reduces 
to (3 + 1) dimensions: N{r, h, 0, (j); Y) Nir, b, 9] Y). The BK equation with 
the impact parameter dependence has been investigated numerically |32[ \'6'6\. 
see also |3J . Here we show some of the results taken from , which used the 
initial distribution in Glauber - Mueller form 



where the impact parameter profile has been chosen to be of the Gaussian type 




iV(°)(r,6,6';y = 0) = I - exp{-r^ S{b)) , 



(29) 



S{b) = l^expi-b^bl). 



(30) 



26 




10.2 Impact parameter dependence 

In Fig. [2ni we show the resulting impact parameter dependence of the BK so- 
lutions calculated for different values of rapidity Y. For small values of b the 
amplitude is large and strong nonlinear effects are evident. On the other hand, 
for large values of b one observes a fast growth of the amplitude is governed by 
the linear part of the equation. One can verify that the increase is exponential in 
rapidity as expected from the BFKL equation. The region in impact parameter 
space, where the amplitude is large, expands with growing rapidity. However, 
perhaps the most striking feature is the fact that the initial profile in impact 
parameter is not preserved even after a small step in rapidity, AY = 0.1. The 
exponential tail of the initial distribution immediately assumes power behaviour 
~ It is interesting to investigate the origin of these power tails. One can 

divide the region of integration into long- and short-range parts by introducing 
a separating cutoff on the dipole size. 



short 



long 



e(ro 



X2 



e(|x2 



ro) 



(Xo - X2)2(xi - X2) 



Y^02 +^^12 ^^01 ^^02 -'*12 J ■ 



In Fig. l27l thc profile in impact parameter space has been decomposed into short 
and long range contributions. One can see that the short range contribution 
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Figure 27: The impact parameter dependence of the amphtude decomposed 
into the short and long range contributions after smaU evolution step in rapidity 
Ay = 0.1. 



dominates the behaviour at small values of b. There the exponential behaviour 
is preserved since we have the factorisation of the initial profile S(b) at small 
values of b. The long range contribution is dominating at large values of impact 
parameter b where it generates the power tail. Thus the ~ 1/6* behaviour 
originates from the integration of the large dipole sizes and is a reflection of the 
asymptotic behaviour of the integral kernel (see discussion in |35j ) . 



10.3 Violation of Proissart bound 

The presence of the power tail in impact parameter has profound consequences 
for the unitarity. Power decrease of the amplitude means that the interaction is 
long range. As already stated it, is a direct consequence of the power-like form 
of the integral kernel 

'2^ ^2 ^2 



X2 



r 



d^xo — r 
54 



^20 ^12 

This type of fast expansion of the interaction system leads to the violation of 
the Froissart bound, as has been first observed in |H5| (compare also a parallel 
discussion in 36 ). It turns out that even though the amplitude is equal or less 
than 1, due to the nonlinearity of the equation, the dipole cross section increases 
fast with the decreasing x which violates Froissart bound 



J (fhN{r,h;Y = Inl/a 
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This happens because the kernel in the BK equation is conformally invariant, 
with no mass scale which would cut off the long range contributions. 



10.4 Dipole size dependence 

The dipole size dependence in the case of the impact parameter dependent BK 
equation is shown in Fig. 1281 where the value for the impact parameter has been 
chosen to correspond to a central collision. For small and moderate values of r 
it has qualitatively the same behaviour as previously, the amplitude vanishes as 
r tends to and extends to lower values of r as rapidity grows. It also saturates 
to 1 for moderate values of r. However, at large values of the dipole size the 
situation is dramatically different. Here the amplitude drops down again. The 
reason is that now there is a dimension in impact parameter which characterizes 
the size of the target. As the dipole grows, at some point it has to completely 
miss the target and amplitude becomes zero again. This is quite different from 
the previous case (without the impact parameter) where the amplitude was 
always saturated since there was an infinite target. It is also interesting to 
study the solutions at different values of b. In Fig. 123 we present the dipole 
size dependence of the amplitude for larger value of h which corresponds to 
a peripheral collision of the dipole with the target. In that case we see that 
the amplitude peaks for values of the dipole size twice the size of the impact 
parameter. This is expected, since it reflects the properties of the integral kernel 
in the BK equation. The solution in the b-dependent case has also several other 
interesting properties. For values of dipole size much smaller than the impact 
parameter r <C 6, the amplitude depends only on a single combined variable 
This dependence on one variable, the anharmonic ratio , is the result of 
the conformal symmetry, see for example 

10.5 Saturation scale with b dependence 

One can also extract the saturation scale from the solution to BK equation in 
(3+1) dimensions. In 32 a following prescription has been used 

{N{r - 1/Q,, 5, 6; Y))e = k, k ^ 0.5 . (31) 

From Fig.|2Slwe see that the above equation possesses two solutions 

<r<RHib,Y). 



Qs{b,Y) 



The lower bound Qs{b,Y) is the impact parameter dependent saturation 
scale, which has been plotted in Fig. 1201 We see that the saturation scale has 
strong dependence on 6, it is largest at small values of impact parameter and 
then decreases for large values of impact parameter. The physical picture is that 
while the impact parameter is increased, one moves from a strongly saturated 
regime to a more dilute one. The tail of the saturation scale is again power like 
^ which is to be expected from the properties of the integral kernel. We 
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Figure 28: Solution of the impact parameter dependent BK equation with 
fixed orientation and position of the dipole for various values of rapidity 
Y = 0.1,5,8,11. 



see that this behaviour is different from the one that could be anticipated from 
initial conditions (compare dashed lines). The saturation scale has following 
behaviour 

Qlib,Y) ~ g{b) exp(a,2A,r), A, ~ 2 , 

where function g(b) is exponentially falling for small values of b and has a power 
like behaviour at large values of impact parameter. 

The second solution R{b, Y) is new compared to the b-independent case, and 
it just reflects the fact the there is an additional scale present, the finite size of 
the target. 

11 Conclusions and outlook. Beyond BK equa- 
tion 

We have described basic properties of the BK equation which is a nonlinear evo- 
lution equation suitable for the description of partonic systems at high density. 
We have shown that the solution to this equation has a property of geometri- 
cal scaling with the characteristic saturation scale. The nice property of this 
equation is the suppression of the infrared diffusion and independence of the 
regularization for the running coupling. The impact parameter dependence of 
this equation leads to the violation of the Froissart bound despite the fact that 
the amplitude is bounded from above. This is a consequence of purely pertur- 
bative approach and the lack of long distance effects such as confinement in the 
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Figure 29: The same as Fig. |2Hlbut for large value of impact parameter 6 = 5. 




Figure 30: Impact parameter profile of the saturation scale for two various 
rapidities Y = 5 and y = 11. Solid lines are the result of the calculation of 
full BK equation. Dashed lines correspond to the exponential behaviour of the 
saturation scale that was expected from initial conditions. 
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BK equation. 

However, the BK equation has been derived by using strong assumptions 
about lack of correlations in the system and it is thus an equation in mean 
field approximation. It is not clear to what extent the BK equation is a good 
approximation to the full Balitsky-JIMWLK equations. In numerical studies of 
the dipole scattering by Salam [37j it has been shown that the fluctuations are 
very important and lead to a very different result as compared with the mean 
field approach. A lot of theoretical effort has been recently devoted to study 
the role of correlations: in [3S| a more quantitative study of the fluctuations 
has been proposed, in [321 a new equation for the generating functional was 
postulated which takes into account correlations in the nuclei; in j^OI a BK 
equation with two absorptive boundaries has been studied; in 01] a numerical 
study of the full JIMWLK equation has been performed for the first time; in 021 
an analytical study of the Balitsky hierarchy restricted to the dipole operators, 
and in 031 El the role of the discreteness of the gluon system and connection to 
the statistical physics have been discussed. We expect that this line of research 
will be continued in the near future and we will be able to understand the 
fascinating and complex theory of strong interactions even better. 
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